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Abstract 


Analytical  development  of  a  set  of  second  order  elliptic  partial 
differential  equations  for  the  generation  of  three-dimensional  curvi¬ 
linear  coordinates  between  two  arbitrary  shaped  bodies  is  presented. 

The  resulting  equations  have  only  two  independent  variables  and  therefore 
require  an  order  of  magnitude  less  working  core  capacity  than  when  equations 
depending  on  all  three  independent  variables  are  considered.  The  method 
also  allows,  in  a  straight  forward  manner,  the  possibility  of  coordinate 
contraction  in  the  desired  regions. 

An  exact  solution  of  the  proposed  equations  for  the  case  of  an  inner 
prolate  ellipsoid  and  an  outer  sphere  with  coordinate  contraction  is  pre¬ 
sented  to  demonstrate  that  by  using  these  equations  it  is  possible  to 
generate  three-dimensional  coordinates  between  analytically  specified 
surfaces  of  simple  forms  by  analytical  means. 

The  fundamental  constraining  equations  which  have  been  adopted  for 
the  generation  of  coordinates  are  ®  and  A^n  *  0,  where  Ag  ls  the 

surface  Beltrami  operator  of  the  second  order. 
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1.  Introduction 


At  present  a  number  of  techniques  are  under  active  development 
for  the  generation  of  three-dimensional  body-oriented  coordinate  systems 
for  use  in  the  numerical  solution  of  the  Navier-Stokes  equations  and 
other  field  equations  where  the  exact  specification  of  the  boundary 
conditions  is  of  prime  importance.  Among  these  efforts  two  easily  dis- 
cernable  groups  can  be  formed,  (i)  algebraic  methods,  and  (il)  the 
elliptic  equations  method.  In  the  first  group  the  grid  points  in  space 
are  obtained  by  some  interpolation  or  blending  functions  scheme  which 
depends  on  the  given  boundary  data.  The  choice  of  the  interpolation 
scheme  or  of  the  blending  functions  is  crucial  in  achieving  a  desired 
order  of  smoothness  and  distribution  of  the  grid  points  in  space.  This 
line  of  effort  has  actively  been  considered  by  Eiseman  [1,2],  Smith  and 
Weigel  [3],  and  Eriksson  [4].  In  the  second  group  of  efforts,  a  set  of 
three  poisson  equations  in  the  curvilinear  coordinates  are  first  inverted 
and  then  solved  for  the  Cartesian  coordinates  under  the  prescribed  values 
at  the  given  boundaries.  Thus  in  essence  all  the  methods  of  the  second 
group  are  a  straight  forward  extension  of  the  work  of  Thompson  et  al  [5] 
in  two  dimensions.  Research  in  this  area  has  been  conducted  by  Mastin 
et  al  [6],  Yu  [7],  Ghia  et  al  [8],  and  Graves  [9]. 

At  this  stage  of  research  it  is  premature  to  compare  the  two  groups 
since  neither  of  them  have  been  fully  investigated  for  their  inherent 
potentials.  However,  based  on  the  success  of  the  differential  equations 
approach  in  two  dimensions,  e.g.  [5],  it  is  desirable  to  further  investi¬ 
gate  the  elliptic  equations  approach  for  the  generation  of  coordinates. 

The  elliptic  equations  approach  presented  in  this  paper  is  different 
from  the  approaches  adopted  in  the  previously  cited  works,  l.e..  References 
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[6]  -  [9] .  The  proposed  method  depends  heavily  on  the  formulae  of 
Gauss  and  on  the  concept  of  principal  curvatures  of  a  surface.  It  has 
been  shown  that  a  fruitful  arrangement  of  the  classical  differential- 
geometric  results  can  yield  a  method  which  is  easily  programmable  on  a 
computing  machine,  and  which  at  any  time  solves  a  two-dimensional 
partial  differential  equation  of  the  form  used  in  Ref.  [5].  In  this 
paper  only  the  theoretical  development  of  the  method  along  with  a  tech¬ 
nique  to  redistribute  the  coordinate  surfaces  near  the  inner  boundary 
surface  has  been  considered.  The  developed  equations  have  been  solved 
for  the  generation  of  three-dimensional  coordinates  between  an  inner 
prolate  ellipsoid  and  an  outer  spherical  surface  in  an  exact  analytic 
form. 
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■  -G?rzig 


2.  Nocation  and  Collection  of  Foraulaa 

In  what  follows,  the  general  coordinates  are  denoted  as  x*  (i  •  1, 
2,3).  However  when  an  expression  has  been  expanded  out  in  full  and  there 
is  no  use  for  an  index  notation,  we  have  set 


c1  -  i  ,  x2  ■  n  ,  x3  -  C  . 


The  derivatives  of  the  position  vector  r  »  (x,y,z)  are  denoted  as 


-1  3X1  *  M 


The  covariant  components  of  the  metric  tensor  are 


8iJ  "  *i  * 


(2.1) 


while  the  contravariant  components  are  given  by 


*1Jg 


(2.2) 


Thus  in  three  dimensions 


g  -  det(g1J) 

*  8llg22833  +  2g12g13g23  ”  ^g23^2gll  “  ^g13*2g22  ”  *g12*  g33 


Writing 


C1  “  g22g33  “  (g23^2 
C2  "  gllg33  "  (g13)2 
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G3  "  811822  "  (g12)2 
G4  “  813823  "  812833 
G5  "  g12g23  "  813822 
G6  "  812813  ~  811823 


we  have 


g11  »  Gx/g  ,  g22  -  G2/g  ,  g33  -  G3/g 
g12  *  G4/g  ,  g13  -  G5/g  ,  g23  «  G6/g 


The  Christoffel  symbols  based  on  the  metric 


8ij 


are 


Fjk  -  gl£[jk,*] 


where 


lik.n  *  rjk  .  r* 


1,3812  38U 

2\xk  ax* 


(2.5) 


(2.6) 


and  repeated  lower  and  upper  indices  imply  summation.  In  the  sequel  we 

have  also  used  the  surface  Christoffel  symbols  which  have  been  denoted  as 

T“  ,  where  the  Greek  indices  range  over  (1,2)  or  (3,1)  or  (2,3). 

PY 

The  coordinate  which  is  held  fixed  to  account  for  the  surface  geometry 


is  denoted  by  a  superscript  in  parentheses.  Thus  the  unit  normal  vector 
on  the  surface  v  ■  const,  is  given  by 


(2.7) 


?<V)  “  <H0  *  X  Hgl 


where 


v  *  1  :  a  ■  2  ,  $  ■  3  (surface  x1  *  const.) 


v-2:a»3,3«*l  (surface  x2  *  const.) 


(2.8) 


v«3:a»l,3«2  (surface  x3  *  const.) 


The  rectangular  Cartesian  components  of  n^  are  denoted  as 


n(v)  *  (X(V)  ,  Y(v)  ,  Z(v)) 


(2.9) 


The  coefficients  of  the  second  fundamental  form  are  denoted  by 
S(v),  T(v)  and  U(v)  defined  as 


S(v)  =  n(v) 


T 


(v) 


U(v>  =  n(v) 


r  (no  sum  on  a) 
~aa 


(v) 

n  *  r  „ 

-aB 


r^g  (no  sum  on  B) 


(2.10) 


where  (v,a,6)  are  In  the  cyclic  permutations  of  (1,2,3),  in  this  order. 

The  partial  derivatives  of  the  second  order  are  expressible  in  terms 
of  the  first  order  as 


r  -  r 
-ij  ij'k 


(2.11) 


For  a  surface  on  which  one  of  the  coordinates  is  fixed,  the  Gauss' 
equations  are 


v'3-J 
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,  tY  _  +  g(v)n(v) 

*aa  aa'y 


-  TY  r  +  x^n^ 
~a3  ae~Y 


(2.12) 


s66  -  ♦  u«.w 


Where  (v,a,|3)  are  in  the  permutational  sequences  of  (1,2,3)  as  shown  in 
(2.8),  and  the  repeated  index  y  implies  summation  on  the  two  indices  of 
a  surface. 

The  sum  of  the  principal  curvatures  of  the  surface  v  =  const,  is,  [10], 


n(v)  Or,  „  Q(v) \  ir 

V  -  2gaeT  +  8ees  )/Gv 


(2.13) 


where  in  writing  equation  (2.13)  for  a  particular  value  of  v,  use  must  be 
made  of  Eqs.  (2.8)  and  (2.10).  We  now  introduce  two  second  order  surface 
differential  operators  by  using  (2.8),  which  for  v  ■  const,  are 


*  "  8638aa  ~  28aB3a8  +  8aa388 


(2.14) 


4V)  5  — [3„{~:(8RR3„  - 
2  a  pp  a  aP  p 

v  v 


+  -  «.8V» 

V 


(2.15) 


As  is  well  known,  the  operator  Aj  is  Beltrami  differential  operator 
of  the  second  order  [10]. 

The  three  space  Christoff el  symbols  which  have  been  referred  to  in 
the  next  section  are  given  below. 
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3.  Formulation  of  the  Problem 


The  principal  idea  of  the  method  to  be  presented  is  to  generate  a 
series  of  surfaces  on  each  of  vhich  a  certain  a 'priori  chosen  variable 
or  coordinate  is  kept  fixed.  Each  surface  to  be  generated  starts  from 
a  given  curve  of  the  inner  body  and  ends  on  the  corresponding  curve  of 
the  outer  boundary,  cf.  Fig.  1.  A  routine,  preferably  a  spline  fit, 
can  then  be  used  to  join  the  successive  generated  surfaces  so  as  to  have 
a  smooth  three-dimensional  computational  net  for  solving  other  physical 
field  equations. 

To  illustrate  the  method,  we  take  v  =  3  or  x3  =  t;  =  constant  on  each 
surface  to  be  generated.  Thus  a  =  1  and  3  ■  2,  viz.,  a  and  B  respectively 
correspond  to  the  coordinates  x1  =  £  and  x2  =  n-  For  the  sake  of  brevity 
of  notation  we  will  not  use  the  superscript  (3)  unless  it  becomes  neces¬ 
sary.  Thus  from  (2.10) 

S  =  n  •  r  ,  T  =  n  •  r  ,  U  =  n  •  r  (3.1) 

-  -  ~5n  -  -nn 


where 


n  *»  iX  +  jY  +  kZ 


Equations  (2.12)  are 


hi  ■  TUEV  +  55 


r.  ”  Tl0r  +  Tn 

-en  12~y 


r  -  Tl~r  +  Un 
-nn  22-y 


(3.2) 


(3.3) 


(3.4) 


(3.5) 
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From  (2.14)  and  (2.15)  the  operators  and  are 


D  =  go9  9  ~  2g. 


32 


+ 


a2 


*22  .  _2  612  acan  Bn  ,  2 

at  an 


(3.6) 


.  .  l  ra_r  i  ,  a  a  w 

L2  -  t^-Ia?{/5-(g22  35  "  812  3d)} 


,  9  f  1  /  9  9  \  ii 

+  3n* /g-  811  3n  "  812  ac)}] 


(3.7) 


We 


now  multiply  equations  (3.3)  -  (3.5)  respectively  by  g22>  ~^g12’ 


g  adding  and  using  equations  (2.13)  and  (2.15)  to  have 


Dr  +  G3(r^A2C  +  r^n)  *  G3?^ki  +  k2^ 


(3.8) 


where 


A2C  “  G3(2g12T12  "  g22Tll  "  811T22) 


A2n  =  G3(2g12T12  ”  822T11  ”  811T22* 


(3.9) 


G3  “  gllg22  ‘  (g12^2 


To  obtain  an  expression  for  k^  +  k2  consider  equation  (2.11)  and 
utilize  the  property  that  n  is  orthogonal  to  r^  and  r^,  so  that 


!  '  ‘  r!l(!  '  -c> 


5  '  %  '  rl2(5  ' 
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5  *  Hnn  “  r22(?  *  *c> 


(3.10) 


where  all  Che  derivatives  with  respect  to  5  are  evaluated  at  5  ■ 
constant.  Multiplying  Eq.  (3.8)  scalarly  by  n  and  using  (3.10),  we 
get 


G3(kx  +  k2)  ■  (n 


^)(S11‘22 


^1?  ”  ^®12^12  + 


*  r3  ) 
»221ir 


(3.11) 


We  now  propose  the  following  deterministic  problem;  Let  £  and  n  be 
the  surface  coordinates  on  the  surface  ?  *  constant,  subject  to  the 
constraints 


A2£  -  0 


(3.12) 


*2n  *  0 

Then  the  Cartesian  coordinates  x,y,z  of  the  surface  satisfy  the  differential 
equations 


Dr  «  G^k^  +  k2^- 


(3.13) 


The  three  scalar  differential  equations  for  the  generation  of  the  Cartesian 
coordinates  are  then 

g22x££  "  2g12x£n  +  8llXnn  “  m  (3.14) 
g22y££  "  2g12y£n  +  8llynn  “  YR  (3-15) 
g22z«  -  2gi2z£n  +  8llznn  “  ZR  (3'16) 


where 
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(3.17) 


r  -  (Xx5  +  Yy;  +  zV(gnr22  -  2g12r|2  +  g22r[1) 

and 

X  .  (y^  - 

Y  ■  <Ve  ■  Vn1'^ 
z  ■  <v„  -  Ve)/,/E3 

Equations  (3.14)  -  (3.16)  form  a  quasilinear  system  of  partial 
differential  equations  in  which  the  components  of  r^  are  assumed  to  be 
known.  Since  the  values  of  x,y,z  are  known  on  the  basic  inner  and  outer 
boundaries  (denoted  at  B  and  °°  respectively  in  Fig.  1),  a  suitable  way 
of  prescribing  r^  can  be  to  take 

r^  -  fi(n)(r^)g  +  f2(n)(r^)<>>  (3.19) 

where  f^(n)  and  f2(n)  are  suitable  weights  having  the  properties 

fl(nB)  *  l  ,  f2(nB)  -  0 
fl(nJ  -  0  ,  f2(n„)  “  l 

For  exposing  the  essential  nonlinear  terms  in  the  factor  R  we  refer 
to  Eqs.  (2.16)  -  (2.18)  in  which  the  r^  terms  have  been  collected 
separately. 

Referring  to  Figure  2,  we  now  solve  Eqs.  (3.14)  -  (3.16)  for  each 
£  ■  const.,  on  a  rectangular  plane  by  prescribing  the  values  of  x,  y  and 
z  on  the  lower  side  (C^)  and  upper  side  (C2)  which  represent  the  curves 


(3.18) 
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on  B  and  «  respectively.  The  sides  Cj  and  are  the  cut  lines  on  which 
periodic  boundary  conditions  are  to  be  imposed.  The  preceding  analysis 
thus  completes  the  formulation  of  the  problem. 


4.  Coordinate  Transformation  (Contraction) 


For  the  purpose  of  generating  coordinates  between  the  space  of  the 
inner  and  outer  boundary  which  can  be  distributed  in  a  desired  manner, 
we  consider  a  coordinate  transformation  from  5  -*•  x  and  n  -*•  o.  Let 


i  -  €<X)  +  50 


n  •  n(o)  +  nB 


(4.1) 


then 


«  -  «o  at  X  *  xo  •  5(xo>  -  0 


n  -  nB  at  o  -  oB  ,  n(oB)  -  0 


Writing 


A(x)  *  ®<°>  -S 


and  denoting  the  transformed  metric  tensor  as  g ,  we  have 


(4.2) 


811  "  gllA2  *  8ii  "  xv  ♦  vl  ♦  «? 


'll  XXX 


gi2  ■  *i2/ex  •  ■  x~x-  ♦  +  *„* 


IT-**  J  J  'V-** 

12  x  0  X  o  X  o 


g22  -  g22/e2  .  8,5  -  *2  ♦  X2  ♦  «? 


‘22  -o  T  7o  T  o 


G3  "  G3/02x2  ’  G3  ”  *11*22  *  **12)4 


(4.3) 


X-X.Y-Y.Z-Z 
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Further  noting  that 
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S.  An  Analytical  Example  of  Coordinate  Generation 


In  this  section  we  shall  consider  the  problem  of  coordinate  genera¬ 
tion  between  a  prolate  ellipsoid  and  a  sphere  with  coordinate  contraction 
near  the  inner  surface.  This  problem  yields  an  exact  solution  of  the 
equations  (4.5)  -  (4.7). 

Let  n  »  n„  and  n  ■  n  be  the  inner  prolate  ellipsoid  and  the  outer 
B  • 

sphere  respectively.  The  coordinates  which  vary  on  these  two  surfaces 
are  t  and  4 .  We  now  envisage  a  net  of  lines  t  -  const,  and  4  ■  const, 
on  these  two  surfaces.  A  curve  on  the  inner  surface  designated  as 
4  -  18 

x  -  coshn0cos4 
a  0 

y  *  sinhn_sin4  cost 
o  o 

z  •  sinhnsint  sint 
B  o 

Similarly,  the  curve  corresponding  to  4 

x  ■  e  cost 

o 


y  ■  e  sincocost 

z  -  e  sln(  sint 
o 

Based  on  the  forms  of  the  functions 
assume  the  following  forms  of  x,y,z  for  the  surface  4  *  4q: 


(5.1) 


4  on  the  outer  surface  is 
o 


(5.2) 


*  -  f(o)cosC0 

y  -  ♦(o)atn?oco8t 

z  -  4(0) sin;  ain£ 
o 

The  boundary  conditions  for  f  and  4  are 

f(oB)  -  coshnB 

n. 

f(a)  -  e 

OD 

♦(oB)  »  ainhnB 

♦  (°»)  *  * 


(5.3) 


(5.4) 


Calculating  the  various  derivatives,  metric  coefficients,  and  all 
other  data  needed  in  the  equations  (4.5)  -  (4.7),  we  get  on  substitution 


an  equation  which  has  sin2Co  and  cos2( 
of  8inzc  and  cos2c  ,  we  obtain 


Equating  to  zero  the  coefficients 

(5.5) 

(5.6) 


where  a  prime  denotes  differentiation  with  respect  to  o.  On  direct  inte¬ 
gration  of  Eqs.  (5.5)  and  (5.6)  under  the  boundary  conditions  (5.4), 

we  get 
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where 


f (o)  -  AeBn(0)  +  C 
♦(a)  -  DeBn(o) 


A 


(e  -  coshnB)slnhnB 
e  -  slnhnB 


B  -  ln[- 


slnhn, 


.jl/(n.  -  nB) 


n» 

e  (co8hnB  -  sinhrig) 
- 

e  -  Binhn_ 

D 


D  «  sinhrig 

As  «n  application  we  nay  take  [11 ] 

C(x)  •  *X 

n(o)  -  b(o  -  oB)Ka 


where  a  and  b  are  constants.  Since  at  n  , 

n(oJ  -  n„  -  nB 


hence 


n(o) 


(n.  -  nB)(o  -  oB)_(o  .  0J 

<«.  -  V 


(5 

(5 

(5 

(5 

(5 

(5 
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By  taking  a  value  of  K  slightly  greater  than  one  (K  ■  1.05  or  1.1) ,  we 
can  have  sufficient  contraction  of  coordinates  near  the  inner  surface. 

For  the  chosen  problem,  since  the  dependence  on  c  is  simple,  we 
find  that  the  coordinates  between  a  prolate  ellipsoid  and  a  sphere  are 

x  -  [AeBn(a)  +  C]cosC 

y  ■  DeBn^a^ sinccosC 

z  *  DeBn^a^sin?sinC 

where  A,  B,  C,  and  D  are  given  in  equation  (5.9). 
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6.  Conclusions 


A  new  method  for  the  generation  of  three-dimensional  coordinates 
between  two  arbitrary  shaped  bodies  has  been  presented.  The  method  is 
based  on  some  simple  differential-geometric  concepts  such  as  the  equations 
of  Gauss  and  the  expressions  for  the  principal  curvatures  of  a  surface. 

The  simplicity  of  the  method  lies  in  solving,  at  one  time,  only  three 
partial  differential  equations  of  the  two-dimensional  type.  This  aspect 
is  bound  to  reduce  the  working  core  requirements  for  a  given  problem  on 
a  computing  machine.  Finally  the  method  allows,  in  a  very  direct  fashion, 
the  possibility  of  coordinate  redistribution  in  the  desired  regions  (cf. 
Eqs.  (A. 5)  -  (4.7)). 

An  analytic  solution  of  the  proposed  equations  for  the  case  of  an 
inner  prolate  ellipsoid  and  an  outer  sphere  has  been  presented.  This 
example  shows  that  one  can  generate  coordinates  between  two  analytically 
specified  surfaces  of  simple  forms  by  exact  solutions  of  the  proposed 
equations. 

In  this  paper  the  fundamental  equations  which  form  a  set  of  con¬ 
straints  for  the  generation  of  coordinates  in  the  surface  are 

&2£  “  0 

a2n  ■  o  , 

where  A2  is  the  surface  Beltrami  operator  of  the  second  order.  It  must 
be  noted  that  A2  is  neither  a  Laplace  operator  in  the  Cartesian  plane 
(x,y),  nor  in  the  Cartesian  space  (x,y,z).  However,  in  the  case  of  a 
Cartesian  plane  (x,y),  when  there  is  no  dependence  on  z,  reduces  to 
the  Laplace  operator  V2. 
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